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We show that oscillating flow of a simple fluid in both the Newtonian and the non-Newtonian 
regime can be described by a universal function of a single dimensionless scaling parameter ojt, 
where a; is the oscillation (angular) frequency and r is the fluid relaxation-time; geometry and linear 
dimension bear no effect on the flow. Experimental energy dissipation data of mechanical resonators 
in a rarefied gas follow this universality closely in a broad linear dimension (10~^ m< L < 10~^ m) 



and frequency (10^ Hz < cj/27r < lO"" 
flows of simple and complex fluids. 



Hz) range. Our results suggest a deep connection between 



The law of similarity P, [if is the most basic conse- 
quence of the Navier- Stokes equations of fluid dynamics: 
two geometrically similar flows with different velocities 
and viscosities but equal Reynolds numbers can be made 
identical by rescaling the measurement units. Similarity 
turns into universality [sl] if additional characteristics of 
the flow - such as the constitutive stress-strain relations, 
the flow geometry, and other dynamical aspects - can all 
be incorporated into the scaling. For steady flows past 
solid bodies, for example, the law of similarity can be 
formulated [2] in terms of the velocity field u: 



u= Uf{-,Re). 



(1) 



Here, LT" is a characteristic velocity of the flow, r rep- 
resents the position vector, L is a dynamically relevant 
linear dimension, and Re = UL/u is the Reynolds num- 
ber with kinematic viscosity u. The dimensionless scaling 
function /, which reflects subtle features of the flow, is 
found from the Navier- Stokes equations. Defining the 
mean free-path and relaxation time as A and r, respec- 
tively, one can express u in terms of these microscopic 
parameters as u ~ ^'^/^ ^ ^^s^ where Cg is speed of 
sound. 

The equations of fluid dynamics can be formulated in 
the most general form as 



duj 
dt 



(2) 



with i,j = x^y^z. The stress tensor cjij can be ex- 
panded in powers of the Knudsen number (Kn = X/L) 
and the Weissenberg number (Wi = r/T) 0, & Q* 
aij = cr-j^ + cr-^^^ + .... Newtonian fluid dynamics, where 
^ij ~ f (f^ ^ with ji = pv^ corresponds to the 

first term of the expansion. To see this, simply consider 
a steady flow past a bluff body of linear dimension L: 
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here. Ma = U / Cg is the Mach number. Thus, the New- 
tonian approximation is only valid when Kn <C 1. With 
A ^ CgT^ it is easy to see that Kn = x ^ "^T ^ ^ 
and ^ ~ ^ 



^ Ma2 • Consequently, Kn <C 1 and Wi <C 1 
correspond to the same physical limit for fixed Ma. This 
statement can also be extended to an unsteady flow past 
a bluff body, where the shed vortices may be thought 
to introduce independent time (Wi) and length scales 
(Kn). Yet, the characteristic oscillation frequency for 
vortex shedding is given hy Vt = ^^^^^ , where St = is 
the Strouhal number. It is well-known that 0.1 < St < 1 
over a wide range of Re depending on geometry. Thus, 
Vtr = ^ ~ T — , and Wi and Kn are Hnked. It is then 
straightforward to conclude that Newtonian fluid dynam- 
ics and emerging similarity relations, e.g. Eq. ([T]), should 
be valid for slowly -varying large-scale flows: Kn <C 1 and 
Wi<Cl. 

Recent advances in micro- and nanotechnology 
[Sl, Isl, biofluid mechanics [l^, rheology [17], and so on 
have resulted in flows at unusual time (Wi > 1) and 
length (Kn > 1) scales, where Newtonian approximation 
breaks down. Many interesting phenomena, such as elas- 
tic turbulence [12], structural relaxation of soft matter 
and enhanced heat transfer in nanoparticle-seeded 
fluids [14], have been observed in this range of param- 
eters. Thus, a universal description of Newtonian and 
non-Newtonian flow regimes is of great importance for 
both fundamental physics and engineering. 

The similarity law given in Eq. ([T]) is valid in the 
Newtonian regime only. For this regime, since both Wi 
and Kn are small. Re is the only relevant dimensionless 
parameter. In the no n- Newtonian regime, the situation 
is different. For a similarity law valid in both Newtonian 
and non-Newtonian regimes, the scaling relation in Eq. 
([1]) could be modified as 



(3) 



This is a prescribed relation including all relevant dimen- 
sionless parameters. In addition to the above-mentioned 
parameters, a dynamic length scale 5 (to be clarified be- 
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FIG. 1: (color online). Flow geometries and dynamical parameters, (a) Laminar steady flow over a semi-infinite flat plate 
with a viscous boundary layer of thickness 6{x) and a free stream velocity U. (b) Unsteady flow generated by an infinite plate 
oscillating at (angular) frequency oo with peak velocity U in the Newtonian limit {cor <^ 1). (c) The same problem in the 
non-Newtonian limit (ur > 1). A wavelength 6^ and a penetration depth 6- emerge instead of the viscous boundary layer 
thickness S. 



low) may enter the scaling function /. The final form of 
the scaling relation is determined by the physical nature 
of the flow. In order to assess the relative importance of 
various dimensionless parameters and reduce the scaling 
relation in Eq. (|3|), we investigate the second order cor- 
rection to the stress tensor derived from the Boltzmann 
equation 0, H, Q: 



^2 dui duj 



(4) 



Considering first the shed vortices behind a body, we 
readily estimate that a^^^ ^ pX^U'^ / L'^ . Thus, a^^^ + 

^(2) ^ p^U ^pX%^ pj,U(l + ^ A) ^ p^^(l + ry 

this simple case, a appears as an expansion in powers of 
Kn or, equivalent ly, Wi. 

In another classic problem - laminar steady flow over 
a semi-infinite flat plate illustrated in Fig. [TJa) - the re- 
lationship between Kn and Wi can still be stated, albeit 
with important differences from above. The key dynamic 
feature of this flow is the viscosity-dominated bound- 
ary layer of thickness 5 = y^^- Here, U is the free- 
stream velocity outside the boundary layer. The New- 
tonian approximation results in = pi^^^ ^ pCgX^. 
With the symmetries in the problem, the contributions 
to the stress tensor of the kind given in Eq. (|4|) disap- 
resulting in a second order correction [4] axy = 
Thus, the expansion becomes 



pear. 



pA%(V-tx) 



Ul 



(1 + ^2). Note that 



the second order term can also be expressed as Wi since 
x/U = T. This example establishes that the relevant 
linear dimension S is no longer a geometric dimension of 
the body but a dynamic characteristic of the flow and 
Kn^ = X/S (X VWi emerges as the scaling parameter. 

Even more unexpected conclusions emerge for the un- 
steady flow generated by an infinite plate oscillating at 
(angular) frequency u with peak velocity and amplitude 
U and A, respectively (Fig. [T](b) and (c)) [15]. In the 
Newtonian limit cc;r <C 1, the Stokes boundary layer 



thickness, 5 = y 77, is the only length scale in the prob- 
lem. In this limit, Kns ex VWi holds as above, since 



Kn^ = A ~ . / ~ ^/uJr. Due to its geometric sim- 
plicity, this problem can be solved in the entire dimen- 
sionless frequency range < cjr < 00 by a summation of 
the Chapman-Enskog expansion of kinetic theory [l^ 
and non-perturbatively [17]. The analytic solution for 
the velocity field is obtained as 



u.iy) = Ue-y/'- cosM - y/5+), 



(5) 



with two new length scales, a wavelength 5+ and a pen- 
etration depth S-: 



f = (l+c.V)V^ 



COS 



sm 



tan 



± 



tan 



(6) 



In the limit cjr 1, disappears from the problem and 
S- becomes the relevant length scale. However, as uor 
00, the first Knudsen number saturates: Kns_ = y > 

^, indicating that in this limit Kns_ cannot appear as 
an expansion parameter in the stress-strain relation or 
as a scaling parameter in Eq. ([3]). The second Knudsen 
number, Kn^^^ = cor, becomes the exact same 

expansion parameter derived from the Chapman-Enskog 
expansion applied to the Boltzmann-BGK equation [l^. 
We conclude that the only relevant scaling parameter for 
the oscillating plate problem valid in both Newtonian and 
non-Newtonian regimes is Wi = cor, re-emphasizing that 
Wi does not contain any information about the linear 
dimensions of the oscillating body. Moreover, as long 
as A <C L and Re ~ 0, the flow over the oscillating 
body remains tangential along the surface, following the 
natural curvatures. 

In order to provide experimental support for the pre- 
dicted universality, we studied energy dissipation in flows 
generated by oscillating solid surfaces. For a large plate. 
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FIG. 2: (color online). 1/Q/ in quartz crystal resonators at 5 MHz, 14.3 MHz and 32.7 MHz as a function of pressure. The 
inset is an illustration of the shear mode. The ujt ~ 1 transition from non-Newtonian (yellow) to Newtonian (blue) flow takes 
place in the white regions at p ~ 60, 200 and 400 Torr. The dotted lines are asymptotes proportional to p. The solid lines are 
fits to Eq. (0 using fitting factors 0.8, 0.6 and 0.9 and S/m = 2.754 m^kg, 9.827 m^kg and 6.296 m^kg (left to right). 



the average energy dissipated per unit time can be ob- 



tained asE = ^f{ujT) 



2 J • j-y \^I^P^y considering the shear 
stress on the plate. Here, S is the surface area and fioor) 
is the scaling function found as Q 



1 



(l+CjV2)3/4 



(1 + uot) cos 



-(1 — uot) sin 



tan ^ UJT 



tan UJT 



■ (7) 



For a mechanical resonator with resonance frequency 
a;/27r, the dissipation can be translated into a fluidic 
quality factor Qf from the relation 



1 

Q~f 



E 




(8) 



where Est = ^^^^ is the energy stored in the resonator 
with mode mass m. In a simple ideal gas, such as nitro- 
gen at pressure p, rcxl/p. can then be expressed 
as a function of p and follows two different asymptotes: 
l/Qf oc p at low p (non- Newtonian) and l/Qf cx p^/^ at 
high p (Newtonian). The transition between the asymp- 
totes takes place at cjr ~ 1, and shifts to higher p as a; 
is increased since r ex 1/p. 

In our experiments, we measured the Q factor of quartz 
crystals, microcantilevers and nanomechanical beams in 
dry nitrogen as a function of p using electrical [18] and 
optical techniques [8|] . We excited the resonators at very 
small oscillation amplitudes {A <C L) around their res- 
onances. The energy losses arising both from the fluid 
and the resonator itself (coupled to the measurement 
circuit) determine the overall (loaded) quality factor as 
^/Qi = ^/Qf{p) + 1/Qr. At low p, l/Qf 0, allowing 
a measurement of l/Qr , and subsequently, 1/Qf{p) [8] . 

In Fig. [21 we show the basic aspects of the universality 
on similar sized macroscopic quartz crystals moving in 
fundamental shear mode (inset) at cj/27r = 5, 14.3 and 



32.7 MHz. In each plot, the slope ofl/Qf{p) changes due 
to the transition from non-Newtonian (yellow) to Newto- 
nian (blue) flow. In order to fit the data to in Eq. ([8|), we 
used experimentally measured S/m values [25|] and the 
recently suggested empirical form r ~ "^'^^^^ — (r is in 
s when p is in Torr) [8] . The flow for the small amplitude 
shear mode oscillations of the quartz crystals matches 
the large plate problem to a very good approximation 
[ii, [20]. The only non-ideality may come from the ve- 
locity distribution on the surface of the resonator [2l|, 
which explains the small deviation of the fitting factors 
from unity. 

The universality requires that the characteristics of 
the flow remain size and shape independent. We estab- 
lish this aspect in Fig. [3l^a)-(c) by comparing data on 
resonators, which span a broad range of linear dimen- 
sion L and oscillate in different modes: a macroscopic 
quartz crystal in shear-mode at 14.3 MH; a microcan- 
tilever and a nanomechanical doubly-clamped beam in 
flexural modes at 1.97 MHz and 23 MHz, respectively. 
We take the dynamically relevant linear dimension of the 
flow as L ^ V^, determined by the surface area. For 
the quartz crystal, Lq ^ 10~^ m, set by the electrode 



diameter 22|, 



23j . For the cantilever, Lc ~ 10~ m and 
for the beam. Lb ^ ^ x 10~^ m. The modes are illus- 
trated in the insets; the flow remains tangential to the 
solid surface along the gentle curvatures. Kn = X/L for 
the devices are shown on the upper axes. Two things 
are noteworthy: all curves look similar, and L or Kn 
appears to have no effect on the flow in our parameter 
range [9]. Fig. [3](d) shows the dissipation data of macro- 
scopic quartz resonators collapsed onto a dimensionless 
plot along with data of smaller resonators from [8]. In 
analyzing the data of flexural resonators, a fitting fac- 
tor of 2.8 was used as opposed to the near unity fitting 
factors used for the shear- mode quartz crystals [26[. 

In our experiments, r is determined by the microscopic 
interaction of gas molecules and a solid surface. By 
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FIG. 3: (color online). Fluidic dissipation across length scales, (a) A micro cantilever with dimensions I x w x t — 125 /xm x 
36 /xm X 3.6 /xm moving in its first flexural harmonic mode in the out-of-plane direction, (b) A quartz crystal of diameter 0.5 
cm and thickness 0.1 mm moving in fundamental shear mode, (c) A nanomechanical doubly-clamped beam with I x w x t — 
17 /im X 500 nm x 280 nm moving in fundamental out-of-plane flexural mode. The upper axes in (a)-(c) show Kn — \/ L with 
L ^ a/S'- The dotted lines are asymptotes proportional to p and the ujt ~ 1 transition takes place at p ~ 20, 200 and 300 Torr. 
All solid lines are obtained from theory, (d) Scaling of fluidic dissipation. The solid line shows the scaling function f{cor) of Eq. 
(O. The dissipation data of macroscopic quartz resonators are collapsed using Eq. (|8]) with the fitting factors and S/m values 
mentioned in Fig. [2l The data points plotted in grey are from smaller flexural resonators of [1]. For the flexural resonators, 
the fitting factors are 2.8; S/m values are calculated from geometry and mode shape 



naively treating the nitrogen as a gas of hard spheres, 
one can obtain r ~ O-isoxio — ^ when p is in 

Torr). This value, however, only reflects interactions be- 
tween gas molecules, as would happen away from surfaces 
in the bulk region of the gas. The observed dependence, 
which is roughly an order of magnitude larger, points to 
the importance of gas-surface interactions. 

Our data provide evidence for a transition from purely 
viscous (Newtonian) to viscoelastic (non-Newtonian) dy- 
namics in oscillating flows of simple gases in simple ge- 
ometries. In our experiments, Re ~ and non-linear 
effects, such as hydrodynamic instabilities and viscoelas- 
tic turbulence, are not present. The observed transition 
is due to the intrinsic dynamical response of the simple 
fluid to high-frequency perturbations. Similar observa- 
tions are commonplace in macroscopic flows of concen- 
trated long-chain polymer solutions, where r can be long 
and, consequently, Wi > 1 due to the relatively slow 
polymer dynamics [III, [HI, [2^ . In rheology, polymers are 
often treated as elastic springs, and viscoelastic behavior 
of polymer solutions is attributed to the direct contribu- 
tion of polymer molecules to the stress tensor. In this 
sense, our work points to a deep dynamical connection 
between oscillating flows of complex and simple fluids 

nil. 
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